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The fluid-dynamical model of a self-gravitating mass of viscous liquid with singular den- 
sity at the center vibrating in fundamental mode is considered in juxtaposition with that 
for Kelvin fundamental mode in a homogeneous heavy mass of incompressible inviscid 
liquid. Particular attention is given to the difference between spectral formulae for the 
frequency and lifetime of /-mode in the singular and homogeneous models. The newly 
obtained results are discussed in the context of theoretical asteroseismology of pre-white 
dwarf stage of red giants and stellar cocoons - spherical gas-dust clouds with dense 
star-forming core at the center. 
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1. Introduction 

Studying highly ideahzed models of vibrating sclf-gravitating fluid masses forms 
one of important issues of theoretical asteroseismology which seeks to explain vari- 
ability in brightness of an observable object as being produced by its vibrations. 
Among central to this domain of astrophysics is Kelvin fluid-mechanical model of 
a heavy mass of inviscid incompressible liquid of uniform densityp undergoing free 
oscillations with nodeless irrotational velocity of fluctuating flow . An outstanding 
importance of this fiducial model, whose extended discussion can be found in Refs. 
[2-7], is that it sets the standard for analytic study of non-radial pulsations of the 
main-sequence stars and serves as an example, according to ChandrasekhaJSl, of "at 
least one problem for which analytic solution can be found" and which illustrates 
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"the type of difRculties one must confront in the other problem", because "in most 
instances, the problems become of such complexity and involve so many parame- 
ters that elementary methods of solution seem impracticable" . It is this model from 
which the very notion of fundamental vibrational mode has come into the theory 
of stellar pulsations^. 

The characteristic peculiarity of a liquid star vibrating in fundamental mode is 
that the velocity of oscillating flow is described by the potential vector field whose 
potential obey the Laplace equation of the form 

V-Mr,i)=0, <5v(r,t) -Vx(r,i), (1) 
V\ir,t)^0, xir,t)^Aer'Peicos9)d{t). (2) 

The last equation exhibits node-free character of the velocity field as a function of 
distance from center to the surface of the star. It is this feature of oscillating flow is 
regarded as the major kinematic signature of fundamental vibration mode. In the 
Kelvin model of homogeneous mass of a heavy incompressible liquid, implying that 
the density is not altered, Sp = 0, the above equations emerge from the continuity 
equation 

5p = -V(/9Jv) = -(5v- V)p- pV • Jv. (3) 

Our purpose here is to explore some peculiarities of the fundamental vibration mode 
in astrophysical object with inhomogeneous distribution of mass. In so doing we 
consider a model of highly inhomogeneous stellar object in which the non-uniform 
density profile has singularity at the star center of the form 

One of the most conspicuous features of this density profile is that the total mass 
of such a star 

/Att 
p{r)dV = —p,R' (5) 

is finite and identical to that for the total mass of the canonical homogenous liquid 
star model. In somewhat different context the model of inhomogeneous liquid star 
with similar density profile has been briefly discussed by ClaytorP. This curious 
feature of the model, which from now on is referred to as the singular star model, is 
interesting in its own right because it permits too analytically tractable solution of 
the eigenfrequency problem for /-mode. Perhaps the most striking feature distin- 
guishing vibrational behavior of inhomogeneous from homogeneous mass is that the 
nodeless non-rotational vibrations of a spherical liquid mass of non- uniform density 
are of substantially compressional character, as it follows from the continuity equa- 
tion. One of the prime purposes of our study here, therefore, is to find out how this 
striking distinction between the density profiles is reflected in the frequency and 
lifetime spectra of fundamental vibration mode. 
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The paper is organized as follows. In section 2, general fluid- mechanical equa- 
tions of nodeless stellar vibrations in fundamental mode are outlined. In section 3, 
the detailed analytic computation of the spectral equations for the frequency and 
lifetime of /-mode in the singular liquid star model is presented followed by their 
comparison with those for the Kelvin fundamental mode in the canonical homoge- 
neous liquid star model which is the subject of Section 4. The newly obtained results 
are highlighted in section 5 and briefly discussed in the context of asteroseismology 
of pre- white dwarf stage of red giants and stellar cocoons - spherical gas-dust clouds 
with dense star-forming core at the center. 



2. General equations of nodeless stellar vibrations in fundamental 
mode 

The state of motion of flowing stellar matter under the action of forces of buoyancy, 

gravity and viscous stress is uniquely described in terms of five dynamical variables, 
to wit, the density p{r,t), three components of the flow velocity v(r, and the 
pressure p{r,t) obeying the coupled equations of fluid-mechanics and Newtonian 
universal gravity 

^ = _ ^ A = ^ + ^ (6) 
dt dxk ' dt dt dxk ' 

dvj _ _dp_ ^ ^dU_ ^ diTik 
dt dxi dxi dxk ' 

dp Tipdp dvk ^ dlnp 

dt p at oxk amp 

V^C/ = -47rGp. (9) 

The equation for density p is the condition of continuity of flowing stellar mat- 
ter. The Navier-Stokes equation for the velocity flow expresses the second law of 
Newtonian dynamics for viscous liquid and equation for pressure is the condition 
of adiabatic behavior of stellar matter; Fi stands for the adiabatic coefficient of 
gaseous pressure in the star. This latter condition means that in the process of mo- 
tion, like vibrations, the time scale of energy exc;hange between infinitesimally close 
elementary volume of liquid is much longer than characteristic period of oscillations, 
see, for example, Refs. [11,12]). The tensor of Newtonian viscous stresses is given 

by 

■Kik = 2r] Vik + (^^ ~ I '^^ '"n ^ik^ '"ik = ^ (VfeVj + \/iVk) ■ (10) 

In computing frequency of fundamental vibration mode, the equilibrium density 
profile p, the pressure Pc at the star center as well as the transport coefficients 
of stellar matter, the shear r] and the bulk viscosities, are regarded as input, in 
advance given, parameters. 
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The potential of self-gravity U{r) and pressure p{r) in motionless, v = 0, state 
of hydrostatic equiHbrium, are the solutions of coupled equations 

V,p(r) = p(r)V, U{r), V^U{r) = -47rGp(r). (11) 

The gravity potential U is determined, in effect, by Poisson equation for internal 
potential Ui and Laplace equation for external one Ue supplemented by the stan- 
dard boundary conditions of the continuity of these potentials and their normal 
derivatives on the star surface 

V^Ui{r)^-ATTGp{r), r < R, V^eW = 0, r > R, (12) 

U.ir)^U.ir) . ^m.-MM . ,13) 

r=R dr dr r=R 

The general solution of equation for pressure is specified by standard boundary 
condition of stress-free surface p{r ~ R) — 0. 

The equations of linear oscillations of stellar matter about stationary state of 
hydrostatic equilibrium of a star with non-uniform equilibriimr density, p = p(r), 
are obtained by applying to ([S])-® the standard procedure of linearization 

p{r) + Sp{r,t), p^p{r)+Sp{r,t), (14) 
v,^v, + Sv,{r,t), K = 0], U ^U{r) + dU{r,t). (15) 

As was stated, in this work we focus on the regime of irrotational vibrations in 
which the velocity field of fluctuating flow subjects to 

V • (5v = 0, V X = 0. (16) 

Given this, the continuity equation takes form 

Sp = -{SvkVk)p- (17) 

The liner fluctuations of the flow velocity are governed by linearized Navier-Stokes 
equation 

pSvi ^ ~V^5p + pV.,5U + SpViU + VkSnik, (18) 

Snik = 277(5wifc, Svik = ^(Vj(5wfe + VkSv,). (19) 
From ([8]) it follows that rate of change in the pressure is controlled by equation 

Sp^-{SvkVk)p. (20) 

Fluctuations in potential of self-gravity inside the star SU = SUi, caused by fluctu- 
ations in density 6p, subject to the Poisson equation 

V^SU ^ -AnGdp. (21) 

The energy balance in the process of oscillations is controlled by equation 

d f pSv"^ 



dt. 



dV = - [((5wfcVfc)5p - p{5vkVk)5U - 5p{5vk'^k)U - 6T:^k5v,k\dV{22) 
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which is obtained after scalar multiplication of ^8|) by 5vi and integration over 
the star volume. To compute the eigenfrequency of vibrations we take advantage of 
the Rayleigh's energy variational principle at the base of which lies the following 
separable representation of fluctuating variables 

5v,{r, t) = a,{v)a{t), SU{r, t) = (^(r)a(t), (23) 

6p{r,t)=p{r)a{t), p(r) = -(afe(r)Vfe) p(r), (24) 

Sp{r,t)^p{r)a{t), p{r) = -(afe(r)Vfe) (25) 

6TT,k{r,t) = 2f]5v,k = 2i]a,k{r)a{t), (26) 

Svik{r,t) ^ aik{r)a{t), aifc(r) = -(Viafc(r), -(-Vfeai(r)). (27) 

Hereafter a.i (r) stands for the time-independent field of instantaneous material dis- 
placements and a(t) for the temporal amplitude of oscillations. The key idea of such 
representation is that it transforms equation of energy balance (|22p into equation 
for a having well-familiar form of equation of damped oscillations 

Ma + Va^-lCa^O, (29) 
M= j p{r)ak flfc dV, ^ j ri{r)a,k a,fe dV , (30) 
IC = [[{akVkp{r)){ak VkUir)) ~ p(r)(afc Vfc)0(r) - (afeVfe)(afeVfe)p(r)]dV(31) 



Here £ is the energy of free, non-dissipative, oscillations and J- is the dissipative 
function of Rayleigh describing their damping by shear viscosity of stellar matter. 
The solution of (|29p is given by 

a{t) = ao exp(— t/r) cosa;(r)t, (32) 

c.'{r)=u;'[l-{u;r)-% c.^ = (33) 

where lu{t) is the frequency of dissipative oscillations damped by viscosity, lo is the 
frequency of free oscillations and r is their lifetime. Thus, to compute the frequency 
and lifetime one need to specify all variables entering integral parameters of inertia 
A4, stiffness /C and viscous friction T>. 

We start with the field of instantaneous displacements afc(r) which is the key 
kinematic characteristics of /-mode. In the star undergoing irrotational oscillations 
in the /-mode, the shape of an arbitrary spherical surface takes the form of harmonic 
spheroids which are described by 

r{t)=r[l + ae{t)P,{C)], C = cos0 (34) 

where Pe{C) is the Legendre polynomial of the multipole order £ specifying the 
overtone number in fundamental vibration mode. In the system with fixed polar 
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axis the potential field of velocity is found from Laplace equation supplemented by 
boundary condition that radial component of velocity on the star surface equals to 
the rate of the surface distortions taking the shape of harmonic spheroids 

VkSvk^O, Svk^S/kX, ^\ = ^ x^Ar'Pi{Oa{t), (35) 
Svrl^^ = Rit), R{t)=R[l + aeit)P,iC)]. (36) 
Taking into account that (5v(r, t) — a{r)a{t) one has 

a(r)=A,V/F,(C), ^^^J^- (37) 

In terms of this field, the linearized equations for the density dp and for the pressure 
Sp reads 

Sp{r,t)^p{r)a{t), p(r) = -(a • V) p(r) = -A,F,(C)^ r^^^^, (38) 

Spir,t)=p{r)a{t), p(r) = -(a • V)p(r) = -A,P,(C)^r^-i^ (39) 

where p{r) and p{r) are the density and pressure of gravitationally equilibrium, 
hydrostatic, configuration. 

To compute variations in the potential of gravity one must consider two equa- 
tions for internal SUi and external 5Ui potentials 

V^dU,{r,t) ^ ~47rGSp{r,t), VV,(r) = -47rGp(r), (40) 
VHUeir,t)^0, V'</)e(r)=0. (41) 

The outlined energy variational method provides a general framework for computing 
the frequency of /-mode in a inhomogeneous Newtonian liquid star with arbitrary 
form of non-uniform density profile. 



3. Fundamental vibration mode in stellar objects with singular 
density at the center 

In what follows we use term singular stellar object or, for brevity, singular star for a 
self-gravitating mass of viscous liquid with the non-uniform and singular in the star 
center density profile given by equation (U whose total mass is identical to that 
for homogenous star model. For such a singular star, the equilibrium potentials and 
fields of universal gravity are given by 



V2[/,(r) = 

^7^(r < R) 
g^{r < R) = 

Ueir > R) 



~47rGp(r), p{r 
207r 



9 



GpsR' 




p{r)dV^—psR\ (42) 



1 - 



3/2' 



GM / r \ 1/2 



i?3 \R 
R^ GM 



^{r > R) 



-VC/e 



GM 



(43) 



(44) 



December 23, 2009 13:31 WSPC/INSTRUCTION FILE MPLA-2009-24- 
40-3257-3274 



Self- gravitating astrophysical mass with singular central density 7 

The hydrostatic pressure obeying the boundary condition of frcc-from-stress surface, 
p{'i')\r=R = 0, reads 

107r„,„_ , „ r\ „ 107r„ 2 2 5 GAf^ 



p{r) = —GpiR{R-T) = P,\l--y P^ = —GpiR' = —-^. (45) 

It is remarkable that the pressure at the star center, where the density has singu- 
larity, is finite and its radial profile is the linear function of distance from the center 

to the star surface. The pressure in the center is defined by equation of state of the 
1 2 

stellar matter . With in advance given equation of state for Pc the rightmost iden- 
tity in (I45|) is considered as definition of the star radius. The internal gravitational 
energy is 

where W is the total gravitational energy of homogeneous star of equivalent mass AI. 
To take into account the compression effect of self-gravity on mechanical property 
of star matter we adopt that radial profile of viscosity coefficient is identical to that 
for the equilibrium pressure profile, that is, of the form 

KO = ^c(l-^) (47) 

where Vc is the shear viscosity in the star center which along with ps and Pc are 
regarded as input parameters of the model. 



3.1. Exact solution of Poisson equation for variations of 
self-gravity potential in singular star 

Having defined the equilibrium profiles of density p{r), the pressure p{r), the shear 
viscosity profile j/(r) and knowing the field of instantaneous displacements a(r) 
we are able to compute the inertia M and viscous friction V. However, in order 
to compute the stiffness /C we must calculate fluctuations in the potential of self- 
gravity 6U{r,t) = (j){r)a{t), that is, to solve Poisson equation for (j){r) with a fairly 
non-trivial right part 

y'MrX) = -fA,Gp,R'/Hr'-'/'P,iO. (48) 
In the spherical polar coordinates we have 

or or sm d OO of) 

- yA,Gp,i?i/2£/-5/2p,(cos0). 



Assuming a solution of the form 



7/ f 7* I 

(j)^{r,e) = ^ Pi{cose) (50) 
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and taking into account that P^(cos^) is the solution of equation 
-^|sin^^^^ = -^(£+l)P,(cos^) 

we obtain 



r^^^ - e(i + i)«(r) = -—Ae Gp, R'/^er'+'/\ 
or'' 3 



(51) 



(52) 



The general solution of equation for , which is finite at the origin, is given by 

Mr,0) = [Cr' + Aep,GR'/'r'-'/'^ P,{cose) (53) 



Outside the star we have 

V'^SUe = ^ V^M^) = o> 
(j,^ = Dr-'^-^Pe{cos0). 

The arbitrary constants C and D are eliminated from boundary conditions 



r=R 



d(f>i _ d^e 
dr dr 



r=R 



which yield 
C = 

Finally, we obtain 



IOtt / r. _ lO'^ 4 



R 



21+1 



(54) 
(55) 



(56) 



(57) 



IOtt £ 
3 2t + l 



Ai psG 



J 2(2^+1)^1 



M + 1 



Pe{cos9). (58) 



It is worth emphasizing that following this line of argument one can get the solutions 
of Poisson equation for a more wide class of inhomogeneous star models undergoing 
node-free spheroidal vibrations with non-rotational field of velocity. 



3.2. Spectral equations for frequency and lifetime of f-mode 

The mass parameter M. is given by 

R 

M = J p{r) [al + al]dV=^p,AjR'/^ J r'^'-^l^ dr 



+1 

/ 



207r ,2 
dC = —A\p,R 



2t+l 



(59) 
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Computation of the viscous friction parameter, with the non-uniform profile of shear 
viscosity (|T7)) . yields 



V 



-ft 



■q{r)r dr 



dC 



+ c 



dC 



if(i-0 



dC 



dC = SnAt £{£ - 1){2£ - 1) / r]{r) r^'^-^dr 



= All A} r/, R^^-\t-l). 
The lifetime, r — 2A^/I?, of ^-pole overtone of /-mode is given by 



3 ''(4^-Hl)(^-l)^ 



(60) 



(61) 



V Ps 

Somewhat lengthy but simple calculation of integral parameter of stiffness /C which 



is presented in Appendix yields 



r - 40^^' ,2^ 2„2£+i e{Ae^ + 2e-i) 
^ - -^^i^Ps^ {2e+i){u + iy 



(62) 



From analytic form of A4 and /C it follows that in the singular star model under 
consideration the lowest overtone of /-mode is of dipole degree, £ = 1. 

The frequency spectrum of fundamental vibration mode in the singular star 
model reads 



L0}{e) = CO, 



2i{2l+\) - 1 



2 47r ^ GM 



2(2^+1) ' 

This last equation can be recast in the following equivalent form 



£> 1. 



Ufit) = UJG 

£ » 1 



£■ 



1 



1/2 



2{2i+l)_ 
tj^(£) ~ cogVI. 



(63) 

(64) 
(65) 



showing asymptotic behavior of the frequency at very high overtones. 

The obtained frequency spectrum of /-mode in singular star model has one and 
the same physical meaning as Kelvin spectral formula for fundamental vibration 
mode in homogeneous liquid star model 

"2£(^- 1)"!^/^ 



UJK{i) = tOG 



[21 +1) 



OJG 



GM 



£>2 



(66) 



but the lowest overtone of Kelvin /-mode is of quadrupole degree, £ — 2. As is 
demonstrated in Fig.l, the most essential differences between the above spectral 



December 23, 2009 13:31 WSPC/INSTRUCTION FILE MPLA-2009-24- 
40-3257-3274 



10 S. I. Bastrukov, H.-K. Chang, E.-H. Wu and I. V. Molodtsova 



6. 



1. 

3^ ' 



0.6 - 



Singular star 
Homogeneous star 



10 



15 



20 



25 



30 



35 



40 



Fig. 1 . Frequency of fundamental modes luj as a function of multipole degree I of global node- 
free irrotational vibrations of singular and homogeneous star normalized to the natural unit of 
frequency luq (upper panel), and the ratio of spectral equations for the frequency of fundamental 
mode in singular model lu j and in the Kelvin homogeneous model cj/f . 



equations are manifested in low-overtone domain (upper panel) and that at large 
£ the frequency spectra of /-mode in singular and homogeneous star models shear 
identical asymptotic behavior (lower panel). In the next section this last spectral 
formula is briefly recovered by the above expounded method with allow for the effect 
of viscous damping of /-mode whose lifetime is computed with non-uniform profile 
of shear viscosity. 



4. Kelvin /-mode in the homogeneous Hquid star 



In the canonical homogeneous star model of uniform density, p — constant the total 
mass has one and the same form as in above singular model, M = (47r/3)pi?'^. The 
internal and external potentials of Newtonian gravitational field are the solutions 
of Poisson equation inside and Laplace equation outside the star 







3GM 






hi 


(Si 




1 - 








2 R 




3 \rJ 



(67) 



2 . 
-VU 



V^U, = 



47r R^ 
U,{t > i?) - — Gp— 
3 r 



GM 



Uir < R) 



GM 



.(r > R) 



GM 



(68) 
(69) 
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In ((68)) . W is the internal gravitational energy W of homogeneous star. The solution 
of equation of hydrostatic for pressure p{r) is given by 



R = 



27rGp2 



(70) 



(71) 



The last expression shows again that the star radius is determined by the central 
pressure Pc related to the densitv o bv equation of state. 




For comparison, in Fig. 2 we plot the fractional pressure and gravity potential 



profiles computed in homogeneous model 



p(r) 



l-(- 
R 



3 GAP 
8^ i?4 



Uo 



and in singular inhomogeneous model 

5 GAP Ui{r) 



Pc y r) ' 



Pr = 



87r i?4 ' Uq 



3 \r) 



5 \R) 



3/2' 



R 
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The governing equations for non-compressional irrotational vibrations of the homo- 
geneous liquid star are 



SiTik = 2r]Svik, Svik = ^(ViJwfc + '^kSvi), rj ^ rjc 

The equation of energy conservation is 

d f pSv"^ 
di 



(74) 
(75) 

(76) 



■dV^- J [iSvkVk)Sp - piSvkVk)6U - 57:,kSv,k] dV. (77) 

On substituting here 

Sv,{r,t)^ai{r)a{t), 6U{r,t) ^ (f>i{r)a{t), 6p{r,t) ^ -{ak\7 k) p{r)a{t) 

we arrive at equation of damped oscillator, Ma + Vd + )Ca = 0, with integral 
parameters defining the frequency and lifetime of the form 



M - 
K. ^ 



pakakdV, V ^2 J ri{r)aik atk dV, 
[{a,Wi) (afeVfe)p(r) + p {ak Vfe)0(r)] dV, 



K_ 



2M 
V ' 



(78) 
(79) 
(80) 



The only unknown quantity is the variations of the gravity potentials obeying the 
Laplace equations 

V^SU, = 0, V'^SUe = (81) 

having the general solutions of the form 

dU^ir <R)^Ce /Pi{C)a{t), 6Ue{r < R) = De r-'^'+^^ Pi{C)a{t). (82) 

The arbitrary constants Ci and Di are eliminated from the standard boundary 
conditions 



UMt)) + 5UMt)) = Ue{r{t)) + SUe{r{t)) 



r=R 



[uMt)) + mrm - mr{t)) + summ 

ar dr 



r=R 



(83) 



(84) 



where r{t) = r[l + a{t) Pi (cos 9)]. Retaining in these equations terms of first order 
in a(t) and putting r ^ R we arrive at coupled algebraic equations for C'e and Di 
whose solution leads to the following final expressions for the time-independent part 
of gravity potential^ 



Air Gp I 



21 + 1 i?^-2 



477 

21+1 



GpR'+^r-^'+^'>Pe{C). 



(85) 
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Computation of integral parameters of the inertia M viscous friction V yields 

A< = 4,A|,«"«^ = J^. (86) 

V = = m„kR'l^^ (87) 

One can consider separately oscillations restored by force associated with gradient 
in fluctuations of pressure and force owing its origin to fluctuations in the potential 
of gravity. In accord with this, the integral parameter of stiffness is written as a 



K, — ICp + ICg 

ICp = - [ (ofeVfe) (a,V,)p(r)dV, ICg = - [ p{a^Vi)(j){r) dV (89) 
Jv Jv 

and, hence, for the squared frequency one has 

where ujp designate the frequency of p-mode and LUg is the frequency of g-mode. For 
ICp and ICg we get 

3 2e + i' (2^ + 1)2 ■ ' 



u'^^l+^l -P = TT> ^l = T7 (90) 




Fig. 3. The squared fractional frequencies of p-mode, g-mode and Kelvin's /-mode as functions 
of multipole degree £, computed in the homogeneous liquid star model. 
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For the frequencies of p-mode and g mode we obtain 

Lol I, ojI^^^ -^iJ^ (92) 



^kU)-^1+^1-^1''-^§^- (93) 



As is illustrated in Fig. 3, is positive, whereas is negative and these sings for p- 
mode and g-mode have one and the same sense as in well-known dispersion relation 
of Jeans uP' — (?JP — AnGp characterizing propagation of longitudinal acous tic wave 
with the velocity of sound Cj in a self-gravitating homogeneous 

fluijm]. 

For the 

lifetime of ^-pole overtone, Tf{£), computed with the non- uniform radial profile of 
shear viscosity rj — •q{r) given by equation (j76p . we obtain 




Fig. 4. The fractional lifetime of £-pole overtones of fundamental vibrational mode in homo- 
geneous and inhomogeneous liquid star models with non-uniform profiles of shear viscosity in 
juxtaposition with Lamb spectral equation for damping time of nodeless irrotational oscillations 
of a spherical mass of viscous liquid of uniform density and constant coefficient of shear viscosity. 



It is worthwhile to compare the computed lifetime spectra with the Lamb spec- 
tral equatiorPl 

"^-"•^(^^ = (2^ + [)(;-i) 

that has been obtained in a similar fashion but assuming that coefhcient of shear 
visco sity has constant value in the entire spherical volume of homogeneous viscous 
liquicilSl. It has been pointed out by JeffreyJ^, however, in the context of geoseis- 
mology that the approximation of uniform viscosity does not allows for the effect of 
self-gravity on mechanical property of matter of astrophysical objects. With this in 
mind, we have supposed that it would be not inconsistent to take radial profile of 
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shear viscosity similar to that for hydrostatic pressure in the star. Also noteworthy 
is that Newtonian law of shear viscosity is equally appropriate for viscous liquid 
and viscoelastic solicJI^. The spectral equations for the time of viscose damping 
of nodeless spheroidal vibrations obtained hear for the first time may be of some 
interest, therefore, for the general seismology of Earth-like planet. In Fig. 4, the 
obtained spectral equations for the lifetime, normalized to t^, as a function of over- 
tone number I are plotted for both singular ([61]) and homogeneous ([94]) models in 
juxtaposition with the Lamb spectral formula (|95p . 



Fig. 5. Artist view of a spherical star-forming object, the stellar cocoon, whose largescale vibra- 
tions in /-mode can be analyzed on the basis considered highly inhomogeneous model. 



5. Summary 

The most striking differences between vibrational behavior of homogenous and in- 
homogeneous liquid star models in the fundamental mode of global nodeless irrota- 
tional pulsations under the combined action of the buoyancy and gravity forces is 
that in the inhomogeneous model the oscillations are of substantially comprcssional 
character, that is, accompanied by fluctuations in density and the lowest overtone of 
/-mode is of dipole degree, as has been first been observed by Smeyers (as pointed 
out in Ref.[8]) and in Ref.[16]. In the homogeneous star model they are characterized 
as non-compressional and the lowest overtone is of quadrupole degree. 

The obtained spectral equations for the frequency and life-time of fundamental 
mode in stellar object with a highly inhomogeneous density profile can be invoked 
in assessing variability of emission from star-forming clouds like stellar cocoon. 
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stellar cocoon 




dusty circumstellar shell 
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pictured in Fig. 5, and too from red giants being on the pre- white dwarf stage as 
being produced by their global nodeless vibrations. These are stellar objects whose 
asteroseismology is adequately modeled by equations of fluid-dynamics, contrary to 
the extremely dense compact objects of finale stage of evolution track, white dwarfs 
and pulsars, the seismic vibrations of which are studied within the framework of solid 
star models relying on equations of sohd mechanics, as reported in Refs. [13,17-19]. 
Understandably, therefore, that mathematical treatment of node- free vibrations of a 
solid star in fundamental mode is different from the considered here fluid-dynamical 
theory of fundamental mode in liquid star models and this difference will be the 
subject our forthcoming paper. 
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Appendix A. Appendices 

The integral parameter of stiffness 

}C = j [{akVk p(r))(afe Vfe f/(r)) - (afeVfe) (afeVfe)p(r) - p(r)(afeVfe).^(r)] dV 
can be conveniently represented in the form /C = Ki + -\- K3, where 

R +1 
IQ^2 r r in^2 



9 

The integral for is given by 
K2 = - j (akVk) {akVk)p{r) ^ = - / 



IOtt^ 
9 



■ ... 5 aeir,e) 8 
+1 +1 



GAR i<-"-'dr en -1)1 P^Qdi + tj{l-C'f" (^) ' < 



407r2 . £2 



21+1 
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In similar fashion, for K3 we obtain 



^3 = - y P{r){ak Vfc)</)(r) «^ = - / 



ar{r,e) 



dr 



+ 



do 



dV 



-507r^ 



GplA^R 



.21-1 



21-Z \ J2 



+ y ( (2ZTT)^ - -uTl) ^ J - ^ ) (^) 
-1 



2^.2 r2£+1 



£ 



9 --'-'--'-^ (2^+ l)(4i'+ 1)' 
The resultant expression for stiffness fC reads 



References 

1. W. Thompson W (Kelvin), Phil. Trans. Roy. Soc. Land. 153, 384 (1863). 

2. H. Lamb, Hydrodynamics (Dover, New York, 1945). 

3. S. Chandrasekharm Hydrodynamic and Hydromagnetic Stability (Clarendon, Oxford, 
1961). 

4. S. Chandrasekhar, Astrophys. J. 139, 664 (1964). 

5. M. L. Azicnman and P. Smeyers, Astrophys and Space Sci. 48, 123 (1970). 

6. S. I. Bastrukov, Phys. Rev. E 53, 1917 (1996). 

7. S. I. Bastrukov, Int. J. Mod. Phys. D 5, 45 (1996). 

8. A. Gautschy, Vistas in Astronomy 41, 95 (1997). 

9. D. D. Clayton, Am. J. Phys. 54, 354 (1986). 

10. C. J. Hansen, S. D. Kawaler and V. Trimble, Stellar Interiors 2nd edn. (Springer- 
Verlag, New York, 2004). 

11. M. J. Thompson, An Introduction to Astrophysical Fluid Mechanics (Imperial College 
Press, World Scientific, 2004). 

12. A. C. Phillips, The Physics of Stars (Wiley, New York, 1994). 

13. S. I. Bastrukov, H.-K. Chang, §. Mi§icu, I. V. Molodtsova and D. V. Podgainy, Int. 
J. Mod. Phys. A 22, 3261 (2007). 

14. H. Jeffreys, The Earth 6th ed. (Cambridge University Press, 1976). 

15. L. D. Landau, E. M. Lifshits, A. M. Kosevich and L. P. Pitaevskii, Theory of Elasticity 
3d edn. (Pergamon, Oxford, 1986). 

16. D. V. Podgainy, S. I. Bastrukov, I. V. Molodtsova and V. V. Papoyan, Astrophys. 
39, 278 (1996). 

17. S. I. Bastrukov, H.-K. Chang, J. Takata, G.-T. Chen and I. V. Molodtsova, Mon. 
Not. R. Astron. Soc. 382, 849 (2007). 

18. S. I. Bastrukov, H.-K. Chang, G.-T. Chen, and I. V. Molodtsova, Mod. Phys. Lett. 
A 23, 477 (2008). 

19. S. I. Bastrukov, G.-T. Chen, H.-K. Chang, I. V. Molodtsova and D. V. Podgainy, 
Astrophys. J., 690, 998 (2009). 



